% 垂直方向 Sobel 计算
[Gx, ~] = imgradientxy(smoothed_image, 'sobel');

% 提取竖线（阈值化）
vertical_edges = abs(Gx) > 30; % 调整阈值以去除噪声

% 保存竖线结果
imwrite(vertical_edges, 'vertical_edges.png');
